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1. INTRODUCTION 

The original Zamolodchikov's C-theorem is related to zero-temperature systems. It 
establishes the existence of a dimensionless function C of the coupling constants with 
monotonic properties along the renormalization group trajectories The assumptions 
presented in the proof are related with the energy-momentum conservation, the 
rotational and translational symmetries, and positivity in a two dimensional quantum 
field theory. The behavior of the C-function reflects the role the quantum fluctuations 
and it is useful in determining the qualitative futures of the theory away from the 
criticality. At the fixed points it takes the value of the central charge of the corresponding 
conformal field theory. Since the basic assumptions underlying the C-theorem are not 
specific to two dimensions only, a considerable interest exists in generalization of the 
Zamolodchikov's result for dimensionalities different from two as well as for nonzero 
temperatures [0] - [^]. Earlier efforts (see p[ and references therein) have been devoted 
to finding a version of the C-theorem valid in four dimensions. There the approach 
was based on a careful investigation of the form of the trace of the energy-momentum 
tensor, written in terms of finite local composite operators. Despite of being able to write 
expressions for the Zamolodchikov's equations for the C-function, similar to the case 
of two dimensions, it turns out that it is not possible to demonstrate the monotonicity 
property. Let us note also the fact that the 3-d analog of central charge is not 
equivalent to the universal number characterizing the size dependence of the free energy 
at the critical point (which is always the case in 2-d conformal field theory). This 
fact indicates that a straightforward generalization of the Zamolodchikov's C-theorem 
is not to be expected for a general d. See also Ref. p, where different approaches to 
the problem have been offered and where it was shown that no direct relations exist 
between the "finite temperature C-theorem" and the Zamolodchikov's C-theorem at 
zero temperature. In the present study an approach, proposed by Netto and Fradkin 
(in some sense a thermodynamic one), for finding a candidate for C-function will be 
considered. In the following dimensionless function is defined 

C (/?, a\d) = -P'^'-^, [f (/?, 9, a\d)-f (oo, g, a\d)] , (1.1) 

where f5 is the inverse temperature [jS = l/T) with the Boltzman constant ks = 1, 
g is a set of dimensionful coupling constants, f {oo, g,a\d) = EQ{g,a\d) is the zero- 
temperature energy density, i.e. the energy of the "infinite" in the inverse temperature 
system, / (/5, g,a\d) is the full free energy density (per unit volume) of the system, and a 
is the characteristic length scale of the lattice. Here n{d) is a positive real number (which 
depends only on the dimensionality d of the system) and v is the characteristic velocity 
(e.g. the velocity of quasiparticles) in the system. The function C (/?, g,a\d) is considered 
to be the d-dimensional nonzero temperature extension of the Zamolodchikov's C- 
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function. It is supposed to be positive, and, in the regions where the quantum 
fluctuations dominate, a monotonically increasing function of the temperature. In 
the numbers n{d) = r((rf + l)/2)C(rf+ l)/7r('^+i)/2 

bosons {C{x) is the Riemann dzeta 
function, r(a;) is the gamma function) and = r((rf+ l)/2)C(rf+ l)(2-2i-'^)/7r('^+^)/2 

for fermions have been suggested. Obviously, the exact choice of n{d) does not effect 
the monotonicity properties of the C-function. Functions analogous to the one defined 
in ( |1.1| ) have been discussed also in Refs. 0] - |@]- 

For a general d the existence of phase transitions in the system, as well as the 
interplay between the quantum and classical fluctuations, makes the analysis of the 
behavior of the C-function difficult from a general point of view. That is why, for any 
d ^ 1 the properties of C have been considered on the examples of concrete models: 
the free massive filed theories (for any d) [Q, the Ising model in a transverse field (for 
= 1) 01 ^^^d the quantum nonlinear sigma model (QNLaM) in the limit N —>■ oo and 
d = 2 0, [^]. (Recently, the value of the C-function at the critical point as a function 
of d , 1 < d < 3, has been calculated in p for that model). 

In the present article we will consider the d-dependence of the monotonicity property 
of the C-function within the framework of an exactly solvable lattice model. We 
will explicitly demonstrate the crucial role that the existence of a quantum (T = 0) 
and/or classical (T 7^ 0) critical points plays for the behavior of C in different regions 
of the phase diagram in the plane temperature - parameter controlling the quantum 
fluctuations. 

The article is organized as follows. In Section II we briefly describe the model and 
in Section III present the basic exact analytical expressions for the free energy of the 
bulk system at nonzero temperature. Section IV contains the analysis of the behavior 
of the C-function in dimensions 1 , 2 and 4 in different regimes of the parametric space 
of the model. Section V presents a finite-size scaling interpretation of the results. The 
article closes with concluding remarks given in Section VI. 

2. The Model 

The model we consider describes a magnetic ordering due to the interaction of quantum 
spins. The Hamiltonian has the following form [Q] 

^ = E - ^ E + f E (2.1) 

where Sg are spin operators at site i, the operators Ve are "conjugated" momenta (i.e. 

= 0, [PiyVi'] = 0, and [Ve,Se>] = i6u>, with h = 1), the coupling constants 
Je,e' = J are between nearest neighbors only, the coupling constant g is introduced so 
as to measure the strength of the quantum fluctuations (below it will be called quantum 
parameter), and, finally, the spherical field /i is introduced to insure the fulfillment of 



the constraint {S}) = 1 . Here (• ■ ■) denotes the standard thermodynamic average taken 
with Ti. 

In the thermodynamic hmit the reduced free energy f^o {P,g\d) = foo {P,g\d) j 
takes the form PI 



A/oo it, \\d) = sup 



t 



X In 



<i> \ (27r) ' 

A 



dqx... / dqd 



:2.2) 



2 sinh 



2t\ 



2^(1 - cosgi 



where we have introduced the notations: A = \Jg/J is the normahzed quantum 
parameter, t = j is the normahzed temperature and = ^ — 2dis the shifted spherical 
field. The supremum is attained at a solution of the mean-spherical constraint that 
reads 
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E 



dqi . 



dqd 



(2.3) 



{2'K)\^t^^J-^ ^-.""> + 2Ef=i(l-cosg,) + &'m2' 

where b = 27rt/A. 

Eqs. (2.3) and ( |2.3| ) provide the basis for studying the critical behavior of the model 
under consideration. 

The critical behavior and some finite size properties of this model have been 



considered in [T^, [[Tj] for 1 < d < 3. Below we present a brief sketch of the derivation 
of the bulk free energy for = 1, 2, 4 at low temperatures. 



3. The free energy at low temperatures 

By using the identities: 

sinh b 1 ~ 



In 



sinh a 2 + vr^m^ ' 



m=—oo 



where ab > 0, a,b are arbitrary real numbers, 

POO 

In (a + 6) = Ina + / exp {—ax) (1 — exp (—6a;)) 
Jo 



dx 

X 



(3.1) 



(3.2) 



where a > 0, a + b > 0, and the Jacobi identity after some algebra at low temperatures 
(j ^ 1) the expression for the free energy (2.3) can be rewritten in the form 



2A/oo (t, \\d) = Aa(0, rf) - (0 + 2d) - As(0, b, d), 



(3.3) 



where 



a(</), d) 



dx 



2-^7? JO x^/^ 



exp [—x(, 



l-(exp(-2a;)/o(2x))'^ + (3.4) 
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b,d) = 2 T — (47rx)~(''+')/' exp (-x0) R ( ^] , (3.5) 
Jo X \xb^ J 

CO 

R{x) = ^ exp (^—xrn^^ , (3.6) 

m=l 

Io{x) is a modified Bessel function, and in ( |3l3| ) is tlie solution of the corresponding 
spherical field equation that follows from ( |3.3| ) by requiring that the first partial 
derivative of the right hand side of ( p.3[ ) with respect to is zero. The above expressions 
are valid for any d. 

In the remainder we will consider the dimensions d = l,d = 2 and d = 4. Then it 
can be shown that: 

a) for d = 1 

,(^.M)=-^V^Er„-'A-,(?^!^) (3.8) 

where 2-^1 is the hypergeometric function and Ki{x) is the MacDonald function (second 
modified Bessel function). 

b) for = 2, and << 1 



where 



a(0,2) 
s(0,6,2l 



a(O,2) + W2(O)0-— 0=^/2 



27r 



— L12 exp 



(3.9) 



Wd(0) 



+ ^Li3 exp 



dqi 



(3.10) 



2(27r)^ 7-^ 



-1/2 



dgd + 2^ (1 - cosg. 



(3.11) 



is a Watson type integral, >V2(0) ~ 0.3214, and Lin{x) is the polylogarithmic function. 

c) for d = 4, and << 1 



a(0, 4) ^ a(0, 4) + >V4(O)0 - K<f)' + ^^0'^', 



(3.12) 
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+^^^Li4 exp - 



^Lis exp 



2Tih 



+ 4;^Li5 exp 



where W4(0) ^ 0.1891 and 

roo 

r= ^ [exp (-2x) Jq {2x)r dx ^ 0.0677. 

JO 



(3.13) 



(3.14) 



Now we have the basic expressions needed to analyze the behavior of the finite- 
temperature C-function as defined by Eq. ( |1.1|) . 



4. The behavior of the C-function 



4.1. The case d = 1 



From Eqs. (|3.3|) , (3.7) and (3.8) it is easy to see that the only nonanalyticity in the 
behavior of the free energy exists at = 0. Then, for < << 1 one obtains, after 



some algebra, that the C-function |15] can be written in the following scaling form 



C(t,A) 



1/4 , 

-Vo exp(-7yo) 



6 



(4.1) 



where the scaling variable is yo = X'^(f)o/t'^, and we have identified v = \fgJ . Here 

00 = 64exp(-47r/A) (4.2) 

is the Haldane gap type behavior of the solution of the corresponding spherical field 
equation for the zero temperature system. Such type of solution is very well known 
from different problems, e. g., one dimensional anharmonic crystal |jl2| and the quantum 
nonlinear 0{N) sigma model in the large N limit (see, e.g., [0, [Q). In deriving ( |4.1| ) 
we have been interested in such a behavior of the nonzero temperature system which 
approaches the corresponding zero-temperature behavior when T — > 0. 

As it is clear from the above expressions, C is a positive and a monotonically 
increasing function of the temperature. 

The behavior of the C-function for = 1 is illustrated in Fig. 1. 



4-2. The case d = 2 

We are interested in the behavior of the C-function around and below the critical 



point only, i.e. < << 1 [|T^, [|TT]. As it is well known the critical point is at 
A = Ac = 1/W2(0) ^ 3.1114 and T = 0. Then, taking into account in Eq. (|0|) that 
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n{2) = C(3)/ {2tc) and v = \^gJ, for the C-functfon we obtain from Eqs. (3.9) and (3.10) 
that C{t, X) = X (x) , where 



X(x) 



C(3) 



x{y-yo) + l{y'/'-yr) 



+ v^Lis (exp (-7^)) + Li3 (exp {-^))] , (4.3) 



with X = TT (1/A — 1/Ac) A/t. Here y = y{x), and yo = yo{x) are solutions of the 
corresponding equations that folfow from (4.3) by requiring that the first partial 
derivative of the r.h.s. of (4.3) with respect to y, and yo, respectively, to be zero. 
These solutions are 

y/y = 2arcsh (]- exp {-2x)] , (4.4) 



2 



and 



— / , A > Ac , . 

The equation (4.3) determines the exact scaling function of C for the case d = 2. 
From the above equations one can see the different behavior of y in the three regions: 
i) renormalized classical, where y tends to zero exponentially fast as a function of x 
{x » 1) ii) quantum critical, where y = 0(1) (for x = 0(1)) and iii) quantum 
disordered, where y diverges as (4x)^ for x << — 1; (y ~ ix'^'^) ^ 5 where x is the 
susceptibility of the system, see [0). The location of these regions is depicted in Fig. 2. 

The behavior of the O-function reflects the existence of these three regions. 

When A < Ac and t ^ 0, from Eqs. (4.3) -(|475|) it follows that 



Cit, A) ^ 1 - ^ exp [-An (1 - A/AJ t''] . (4.6) 

One explicitly observes the exponentially small corrections to the limit value of O = 1 
(at t = 0) that corresponds to massless bosons in d dimensions 0, 0. 

At A = Ac, C{t, A) simplifies, by using the Sachdev's identity Q], and becomes 

0(t,Ac) = | (4.7) 



This universal rational number |]T6[ has been derived for the first time for the quantum 
nonlinear 0{N) sigma model in the limit oo [^]. It demonstrates that at the 

quantum critical point A = Ac the O-function does not depend on the temperature. The 
difference from the corresponding results in Ref. (compare Fig. 3 in 0] with Fig. 3 
in this article) is due to the fact that terms proportional to the difference between y and 
yo in (4.3) have been neglected there. The above is justified when y >> 1 (then y and 
yo are exponentially close to each other). The analysis of the corresponding equation 
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shows that the last happens when x « —1 (i.e. A > Ac) where y ~ yo ~ (4x)^, which 
is the case of the quantum disordered region. In this case it is easy to see that 

C{t, A) ^ i^^^^ exp [An (1 - A/AJ r'] , (4.8) 

i.e. C approaches zero exponentially fast in terms of the scaling parameter x. The 
behavior of the C-function in this case is that one of massive free bosons |^. 

Let us consider now the monotonicity of the C-function. From (4.3) it follows that 

^^ = -^(l/-yo)(l-A/AJt-. (4.9) 

Since y > yo, we conclude that C is a monotonically increasing function of the 
temperature for A > Ac, and monotonically decreasing function for A < Ac. Within 
exponentially small in temperature corrections this result coincides, in fact, with the 
corresponding one for the QNLcxM in the limit ^ oo 0. 

The above results for the behavior of the C-function are illustrated in Fig. 3. 

Finally, it seems worthwhile to mention that, as follows from Eq. (4.3), the C- 
function is monotonically increasing function of the scaling variable x (see Fig. 4) for 
any value oft (A/t >> 1). 



4-3. The case d = A 

For this case, taking into account that in Eq. ( |1 . 1| ) n(4) = 3^(5)/ (27r)^ and v = \fgJ 
for the C-function we obtain from Eqs. ( |3.12| ) and (3.13) that C(t, A) = X(x, A/t), 
where 

xiy ~ yo) + jr (y^ - yl) - (4.10) 



X{x,X/t) - ^^"^ 



3C(5) 



+ T^Al [^/Lis (exp (-yy)) + 3v^i4 (exp (-v^)) + 3Li5 (exp (-^/y))] 
(27r) 



with 



and Ac = 1/W4(0) ~ 5.2882. Here ?/ > 0, and > are solutions of the corresponding 
equations that follow from (4.10) by requiring that the first partial derivative of the 
r.h.s. of (4.10) with respect to y, and y^, respectively, to be zero. This leads to the 
following equation for y 

^ = + 2^ ^^^^ ^'^^^ ^ ^^^^ (- W))] • (4-12) 
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It is easy to see that, for a given t and A, the solution of the above equation, if it exists, 
is unique. For i/q we get 

-{2t)/{rX)x, 



-X 



yo 



0, 



A > Ac 
A < A, 



(4.13) 



One observes that in the most general case the function X, given by Eq. (4.10), could 
not be recast in a scaling form. However, as we will see below, the last is possible in 
some subregions of the A — t plane. We recall that the susceptibility x of the system is 
proportional to y^^ (if y 7^ [1^, [1^), which leads to the conclusion that a nonzero- 
temperature phase transition exists at a given = tc{X), where tc{X) is given by the 
equation 

.1/3 



tc(A) = A 



1 



(4.14) 



C(3) 

(at t = tc{X) one has y = 0, and y = also for t < tc{Xj). As for the d = 2 case 
three principal different regimes exist: i) renormalized classical (where y tends to zero 
exponentially fast as a function of A/t), ii) quantum critical (where y tends to zero 
algebraically as a function of A/t or y = 0(1)), and iii) quantum disordered (Vhere 
y » 1). In order to describe the behavior of the C-function below we analyze these 
three regimes. 

A) Let us first suppose that y <<1. Then Eq. ([4.12|) becomes 

^ C(3) 1 „A_, ^^^^^1 



1 

a: 



(27r) 



—-2?/ In {y/e)-r-y. 
(4vr) 



A 

i' 



Obviously, there are two subregimes a) when the first term in the right hand side 
dominates and b) when the second one dominates. The borderline between them is 
given by 




1 
A 



(2vr) 









= exp 



- {Anfr^ + 1 



In the X — t plane Eq. (|4.16|) determines a line t*{X) that is exponentially close to (as a 
function of A/t ) the line t^ (A) . At t*{X) the solution of Eq. j ^A^ is 



y ~ exp 



(4.17) 



whereas y = at tc{X). We conclude that the renormalized classical regime is observed 
for parameters of the system laying in the X — t plane between tc{X) and t*{X). In 
this regime one could neglect the second term in the r.h.s. of (4.10) which leads to a 
scaling form of the C-function with a scaling variable x, defined in ( ^.11| ). At t*{X) the 
C-function could not be rewritten in a scaling form. In the remainder we will see that 



10 



another scaling variable x can be defined for the region to the right of t*{X). This is 
due to the fact that in this region the first term in the r.h.s. of Eq. ( [4.12| ) is of the 
leading order. Indeed, this is true not only for the case h) but also for the cases B), 
when y = 0(1), and C), when y ^ 1, which cases are to be considered below. 

Before passing to the consideration of case B) let us note that the further inspection 
of Eq. ( |4.15| ) for the case b ) leads to the conclusion that there exists a crossover line 

1/3 



ts{\) = A 



{2Tlf / 1 1 



(4.18) 



,C(3) VA, A. 

between two regimes where x(^) A) ~ t^^ and x{^) A) ~ respectively. This curve is 
symmetric to the curve tc(A) with respect to the hne A = Ac- Let us turn now to the 
case 

B) y = 0(1). Then, since t << 1, Eq. ( |4.12|) becomes extremely simple and, up to 
the leading order coincides with the corresponding equation for the zero-temperature 
system (see Eq. ( |4.13| )). Its solution is 

Since xi^, A) ~ (yt^y and y = 0(1), one concludes that in this regime x(^5 A) ~ t ^. In 
a given sense a formal curve ti(A) in the X — t plane which borders the region in which 



y = 0(1) can be obtained by simply setting y = 1 in Eq. ( 4.19 ). Summarizing the 



results from A) and B) we are lead to the conclusion that the quantum critical regime 
is observed for values of the parameters t and A laying in the X — t plane between the 
curves t*(A) and ti(A). We see also that in this region X{t, A) = X{x), i.e. the scaling 
is restored with the new scaling variable x. We pass now to the case 

C) y » 1. Then one formally receives the same solution as given by Eq. ( |4.19| ) but 
now x{t, A) ~ (1/Ac — 1/A)~^ , i.e. it does not depend on t up to exponentially small 
in X/t corrections. We conclude that the region of parameters in the X — t plane below 
ti(A) determines the quantum disordered region. Again, as in B), the scaling variable of 
C is X. 

The above results are summarized in Fig. 5. 

The existence of the regions of thermodynamic parameters determined above is 
reflected by the corresponding behavior of the C*-function given by Eq. (4.10). 

First, for t < tc(A) (i.e. under the existence of a long-range order in the system), 
since y = y^ = 0, one immediately obtains from (4.10) that = 1 [pO|] . 

Further, for t > tc{X) taking into account Eqs. (|4.13|) and ( |4.12D , it is easy to see 
that 



dCjt, A) _ 27rV A 



{y - yo) 



x + lr{y + yo] 



(4.20) 
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Since ?/ > yo > for the considered region of parameters A and t, the above equation 
leads us to the conclusion that for any A < Ac the C-function is a monotonically 
decreasing function of the temperature. The same is true also for A = Ac and t 7^ 0. 
From the analysis of the solutions of the equations for y and it becomes clear that for 
a fixed A and t > ts(A) the leading-order form for both y and yo is y = yo = —x (if one 
takes into account next-to-the-leading order terms then, of course, y > yo). Setting the 
above expressions for y and yo in the rectangular brackets in Eq. ( |4.2(J| ), we conclude 
that C is a monotonically increasing function of t for any t > ts{X). It is clear that 
somewhere between A = Ac and ts(A) the derivative of the C-function changes its sign, 
i.e. there is a line tst(A) of stationary points dC{t, X)/dt = . One can see that 

.1/3 



tstiX) = A 



{AttY f 1 1 



(4.21) 



C(3) VAc A. 

Since the point A = Ac,t = lies on tst(A) and at it C = 1, we conclude that C = 1 at 
the whole line tst(A). It is clear now that C is a monotonically increasing function of t 
for t < tst(A). It is a monotonically decreasing function of t for t > tst(A) as well as for 
any (small) t if A < Ac. These results are summarized in Fig. 6. 

5. Finite-size scaling interpretation 

It is interesting to interpret the bulk critical behavior of the C-function in the context 
of the finite-size scaling (FSS) theory by introducing a finite "temporal" dimension 
= X/t. Then, taking into account: i) the dimensional crossover rule that connects the 
properties of a given (i-dimensional quantum system and those ones of the corresponding 
{d + l)-dimensional classical one with the mapping — > d + 1, L ^ L^, t X, ii) the 
Privman and Fisher |^ hypothesis for the free energy of a finite classical system when 



the hyperscaling holds (i.e. between the lower di and the upper du critical dimensions 
of the system), one could make the statement that the free energy f^o of a quantum 
system with dimensionality di < d < du should have the form 

[f^ (t, X\d) - f^ (0, X\d)] /(TLr)-' = L-^^+^^Y {L^/i (0, A)) , (5.1) 

where ^ (0, A) is the correlation length of the zero-temperature system, and F is a 
universal function. We remind that in order to have no nonuniversal prefactor in front of 
Y for classical systems one considers f^o = Pfoo, instead of /oo itself. The normalization 
of the free energy in ( |5.1| ) simply follows by our choice of L^. For the model considered 
here the inspection of Eq. (4.3) shows that the hypothesis ( p.l| ) is indeed valid with 
the standard scaling variable Lr/^ (0, A) = x = vr (1/A — 1/Ac) A/t and C (t, A) = X{x) 
= —Y{x)/n{d). It is interesting that despite the lack of hyperscaling at c?/ = 1 the 
C-function again can be written as a function of Lt-/^(0, A) (see Eq. ([4.1D), if one 
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identifies ^ (0, A) = 0q , where 0o is given by Eq. (^4.2| ). The case d = A > du is much 
more interesting due to the lack of hyperscaling. In the most general case the C-function 
could not be even recast in a FSS form (see Eq.(4.10)). The last is possible exponentially 
close (in L^-) to the line of finite-temperature phase transitions tc (A) ,where the modified 
scaling variable is 2x = [Lt-/^ (0, A)]^ (see Eq. ( [4.11| )). The standard scaling variable 
X = Lr/C, (0, A) is restored only for parameters to the right of the curve t* (A) in the 
A — t plane (see Eq. ( |4.19| ) and the comments connected with it). This change of scaling 
variables from x to x is a new point within FSS theory. Normally one observes modified 
FSS (see [^) above du due to the existence of a dangerous irrelevant variables 
in the system. On the other hand, considering the d = 5 dimensional spherical model 
film Barber and Fisher, as early as in 1973 ||2^, stated that the scaling variable should 
be the standard one, i.e. L^-fC, (0, A) in our notations. The above results resolve this 
seeming contradiction: the scaling variable has to be modified very close to the phase 
boundary, but is the standard one a bit away of it. The physical reasoning for that 
difference is the existence in the system of a temperature driven phase transition in 
addition to the quantum one with respect to A at t = 0. To our knowledge all other 
examples considered previously in the literature of modified FSS concerns finite systems 
with no (sharp) phase transition in it. 



6. Concluding remarks 

One generally expects that the C-function increases monotonically when the quantum 
fluctuations "dominate" 0. The real meaning of the term "dominate" turns out to be 
quite subtle, as we have demonstrated in the current article. In fact, we have shown 
that the region where the C-function remains monotonically increasing (as a function 
of temperature) and the quantum critical region do essentially intersect but do not 
coincide (see Fig. 2). This is one of the results of the present work. The question of 
where one should look for and what should be understood as domination of quantum 
fluctuations is, indeed, very intriguing. It is a part of the more general problem of a 
quantitative description of the interplay of the quantum and critical fluctuations. There 



exist different views on that issue. The standard one [21], |22| is based on the "ratio" 
between the correlation length and the length of De Brougle. Another possible approach 
can be based on the behavior of the C-function 0,0. Furthermore, there is an approach 
based on the algebra of critical fluctuation operators, due to Verbeure and Zagrebnov 
23| , where a measure of the "degree of criticality" is introduced in a mathematically 



rigorous way. 

In the present work we investigated the behavior of the C-function for d = 1,2, 4. 
The case d = 1 represents the situation with no phase transition and strong quantum 
fluctuations, d = 2 - the one when a quantum critical point appears at T = 0, and 
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d = 4- when there is a hne of classical critical points ending up with a zero temperature 
(quantum) critical point. In fact, these are the most typical cases on which the attention 
in the literature is focused. 

Case d = 1. As it is to be expected on general grounds, the C-function increases 
monotonically as a function of temperature (see Fig. 1). This reflects the fact that the 
quantum fluctuations are strong enough (as it is clear from Eq. ( [4.2| ) one cannot consider 
A as a small parameter) and the lack of a critical point. The C-function obtained here 
coincides with the C-function of the massive free bosons (for d = 1) with mass a/0o 
, because of the exponentially small difference then between (p and 0o, i.e. one can 
consider as a fixed parameter in (3.7) and in (3.8). The general case (for any d) of 
free massive bosons actually follows from (p.3|) - ( p.6|) by considering there as a fixed 
parameter connected to the mass m of bosons (0 ~ ). For the last case it is trivial 
to check that the corresponding C-function is that one obtained in (see Eq. (3.2) 
there). 

Case d = 2. Fig. 2 shows the phase diagram for our model which coincides with 
the phase diagram of the d = 1 quantum Ising model, as well as with the nonlinear 
0{n) sigma model in the limit n oo, see, e.g., Sachdev ||21|. As a function of the 
temperature C is monotonically increasing for A above Ac, equals 4/5 at A = Ac (and then 
C does not depend on t) and is monotonically decreasing for A below Ac (see Fig. 3). The 
lack of overall monotonicity with respect to the temperature is due to the crossover from 
classical to quantum behavior. It is clear, that one indeed can consider monotonicity of 
C as a measure of the role the corresponding fluctuations are playing in a given region 
of parameters. It is interesting that C changes its monotonicity, in fact, in the middle 
of the quantum critical region. Finally, we note that it is nevertheless possible to find 
a (nontrivial) variable, with respect to which the C-function is monotonic in the whole 
t— A plane (see Fig. 4). This variable is the scaling variable x = vr (1/A — 1/Ac) A/t. 

Case d = 4. The existence of a line of non-zero temperature critical points modifies 
drastically the corresponding picture in comparison with the d = 2 case. Now a line of 
stationary points tst(A) appears (see Fig. 5) which "starts" from (A = Xc,t = 0) and lies 
to the left of ts(A). To the left of tst(A), C is a non-increasing function of the temperature 
(see Fig. 6). For A < Ac and t < tc(A) one has C = 1, whereas within the region 
between tc(A) and t^* (A) the C-function is monotonically decreasing as a function of the 
temperature. To the right of tsf(A) the C-function becomes monotonically increasing 
function of the temperature, being zero at the t = 0, A > Ac line. At the lines tc(A) 
and tst{X) the C-function reaches its maximum value, i.e. it becomes C = 1. Finally, we 
would like to mention that, similar to the case c? = 2, it is possible to find two nontrivial 
parameters such that with respect to both of them the C-function is monotonically 
increasing. Such are, e.g., the parameters x and X/t (see Eq. (4.10) and take into 
account that y > i/q). 
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Comparing the behavior of the C-function for d = 1, d = 2 and ci = 4 we conclude 
that 

a) For d = 1 for any fixed A we have a monotonically increasing with temperature 
C-function. 

b) For d = 2 the above is true only for A > Ac. 

c) For d = 4 the C-function is a monotonically increasing function of t for A > Ac 
and t small enough. The monotonicity of the C-function does not change by increasing 
t only for A = Ac- 

So, the region in the parametric space where C remains monotonically increasing 
with t becomes smaller when d increases. We explicitly see the crucial role the 
dimensionality d and the existence of phase transition, which appears upon increasing d, 
play in the behavior of the C-function as a function of t. Nevertheless, for any d one can 
find nontrivial variable(s), function(s) of the temperature and the parameter controlling 
the quantum fluctuations, in terms of which C is a monotonically increasing function 
of its variable(s). In close vicinity of the quantum critical point the C-function is given 
by a universal scaling function which properties can be interpreted in terms of FSS that 
has to be modified for d = 4. 
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Figure Captions 

Fig. 1. The bold line = 8Aexp (— 27r/A) borders from above the region in the t— 
A plane where the expression for the 1-d C-function, given by Eq. ( |4.1| ), is valid. The 
symbol C t means that C increases in whole that region starting form C = at t = 0. 

Fig. 2. The phase diagram of the model and crossovers for the case ci = 2 as a 
function of t and the normalized quantum parameter A. One distinguishes renormalized 
classical, quantum critical and quantum disordered regions. Long-range order is present 
only at t = for A < Ac. 

Fig. 3. The behavior of the 2-d C-function is illustrated for A = 1.5Ac, A = Ac and 
A = 0.5Ac. 

Fig. 4. The behavior of 2-d C as a function of the scaling parameter x = 
7r(l/A-l/Ac) X/t . 

Fig. 5. The phase diagram of the model and crossovers for the case d = 4. 
Long-range order exists below the line tc- The line tst is the locus of points in the 
thermodynamic space where dC/dt = 0. The other lines denote crossovers between 
different regimes which are described in the text. 

Fig. 6. The monotonic behavior of C as a function of temperature is shown. The 
symbols C t ( C |) mean that C is a monotonically increasing (decreasing) function of 
the temperature. The number in ellipses show the value of C at the corresponding line. 
In whole long-range order region, i.e. below the line t^, C = 1. 
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